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Abstract 



A detailed analysis of the radiative pp annihilation is made in the frame- 
work of a two-step formalism, the pp annihilates into meson channels con- 



taining a vector meson with a subsequent conversion into a photon via the 

vector dominance model (VDM). Both steps are derived from the underlying 

quark model. First, branching ratios for radiative protonium annihilation are 

calculated and compared with data. Then, details of the isospin interference 

are studied for different models of the initial protonium state and also for 
• i— i . 

different kinematical form factors. The isospin interference is shown to be 

a: 

uniquely connected to the pp — nn mixing in the protonium state. Values of 
the interference terms directly deduced from data are consistent with theo- 
retical expectations, indicating a dominant pp component for the 1 Sq and a 
sizable nn component for the 3 Si protonium state. The analysis is extended 
to the pp — ► 7$ transition, where the large observed branching ratio remains 
unexplained in the VDM approach. 

Keywords : nucleon-antinucleon annihilation, vector-meson dominance, radia- 
tive decay 
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I. INTRODUCTION AND MOTIVATION 



Nucleon-antinucleon (NN) annihilation, due to the richness of possible final meson states, 
is considered one of the major testing grounds in the study of hadronic interactions. Both 
quark and baryon exchange models have been applied to NN annihilation data. 
However, the task of extracting information on the dynamics of the NN process is enor- 
mously complicated by the influence of initial and final state interactions. Some of the 
simplest annihilation channels, where the theoretical complexity of the NN annihilation 
process is partially reduced, are radiative two-body decay modes, where final state inter- 
action is negligible. Experimental branching ratios for radiative decay channels in annihi- 
lation from pp atoms were made available by recent measurements of the Crystal Barrel 
collaboration at CERN, performing a systematic study of the reactions pp — > 7X where 
X = 7 , 7T° , i] , u and rj ||. Radiative decays of the pp atom where, in contrast to ordinary 
production of nonstrange mesonic final states, isospin is not conserved, are well suited for 
studying interference effects in the isospin transition amplitudes |5||J. 

The simplest and most natural framework in studying radiative decay modes is the vector 
dominance model (VDM) j7|. In setting up the annihilation mechanism one adopts a two- 
step process where the pp system first annihilates into two mesons, with at least one of the 
mesons being a vector meson (p and u), and where the produced vector meson converts 
into a real photon via the VDM |J. In this case, production rates of radiative decay modes 
can be related to branching ratios of final states containing one or two vector mesons. A 
first analysis || in the framework of VDM was performed by Crystal Barrel, showing that 
the interference in the isospin amplitudes is sizable and almost maximally destructive for 
all channels considered. The phase structure of the interference term is determined by two 
contributions: i) the relative signs of the generic strong transition amplitudes for pp — > Xuo 
or Xp acting in different isospin channels; ii) the presence of the initial state interaction in 
the pp atom, which mixes the pp and nn configurations. Similarly, analogous sources are 
responsible for the isospin interference effects in the strong annihilation reactions pp —>■ KK 
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|I]j8|||. Here, however, definite conclusions concerning the size and sign of the interference 
terms depend strongly on the model used for the annihilation dynamics. 

In the present work we show how the determination of the interference terms in the 
analysis of the radiative decays can be uniquely connected to the isospin mixing effects in 
the pp atomic wave functions. The extraction of the magnitude and sign of the interference 
from the experimental data can in turn be used to investigate the isospin dependence, at 
least in an averaged fashion, of the S-wave NN interaction. We study this point for different 
NN interaction models. 

This paper is organized as follows. In Sec. |I| we develop the formalism for radiative 
decays of protonium. As in Ref. || we adopt a two-step formalism, that is pp annihilation 
into two-meson channels containing a vector meson and its subsequent conversion into a 
photon via the VDM. Both steps are derived consistently from the underlying quark model 
in order to fix the phase structure of the isospin dependent transition amplitudes. We also 
indicate the derivation of the branching ratios for radiative decays of S-wave protonium, 
where the initial state interaction of the atomic pp system is included. Sec. |TJ is devoted to 
the presentation of the results. We first perform a simple analysis to show that theoretically 
predicted branching ratios for radiative decays are consistent with the experimental data. 
We then show that the isospin interference terms present in the expression for the branching 
ratios can be uniquely connected to the pp - nn mixing in the atomic wave function, induced 
by initial state interaction. We quantify the details of this effect for different models of the 
NN interaction and apply the formalism developed in Sec. [lj to extract size and sign of the 
interference from data, which will be shown to be sensitively dependent on the kinematical 
form factors associated with the transition. Furthermore, we comment on the application 
of VDM on the transition pp —>■ 7$, where the corresponding large branching ratio plays a 
central role in the discussion on the apparent violations of the Okubo-Zweig-Iizuka (OZI) 
rule. A summary and conclusions are given in Sec. [IV]. 
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II. FORMALISM FOR RADIATIVE DECAYS OF PROTONIUM 



In describing the radiative decays of protonium we apply the vector dominance model 
||[7]]. We consider the two-step process of Fig. 1, where the primary pp annihilation in 
a strong transition into a two- meson final state, containing the vector mesons p and u, is 
followed by the conversion of the vector meson into a real photon. Here we restrict ourselves 
to orbital angular momentum L=0 for the initial pp state, corresponding to the dominant 
contribution in the liquid hydrogen data of Crystal Barrel Furthermore, we consider 
the transition processes pp — *> jX, where X = 7, tt°, 77, p, u and rj', with X = presently 
excluded. The final state 07 plays a special role in the discussion of the apparent violation of 
the Okubo-Zweig-Iizuka (OZI) rule, where a strong enhancement relative to w~f was observed 
|T0| . Within the current approach the description of the first-step process pp — ► cj(p)0 and 
its phase structure cannot be accomodated due to the special nature of the 0, a dominant ss 
configuration. Later on we will comment on the possibility to explain the enhanced 07 rate 
within the VDM, as suggested in the literature [|ll|], and on the implications of the analysis 
presented here. 

In the two-step process we have to introduce a consistent description for either transition 
in order to identify the source of the interference term. In particular, the relative phase 
structure of the strong transition matrix elements pp — > uM versus pp — > p°M (M = 
7T°, i], p, u and 7/) is a relevant input in determining the sign of the interference. Basic SU(3) 
flavor symmetry arguments [|l^] do not allow to uniquely fix the phase structure, hence 
further considerations concerning spin and orbital angular momentum dynamics in the NN 
annihilation process have to be introduced. Microscopic approaches to NN annihilation 
either resort to quark models (for an overview see Ref. |J) or are based on baryon exchange 
models |2|-44| . Here we choose the quark model approach, which allows to describe both, the 
strong transition of pp into two mesons and the vector meson conversion into a photon. 

For the process pp — > VM where V = p, u and M = it , 77, p, uo and 77' we apply the 
so-called A2 model ffl, depicted in Fig. 2a. In the discussion of annihilation models based on 
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quark degrees of freedom this mechanism was shown to give the best phenomenological de- 
scription in various meson branching ratios [|I|,|1^,|1^] . In a recent work |]15[ we showed that 



the A2 model combined with a corresponding annihilation mechanism into three mesons 
can describe pp cross sections in a quality expected from simple non-relativistic quark mod- 
els. The transition matrix element of pp — > VM in the A2 model including initial state 
interaction is given by: 

T N n { ij)^vm = < V(ji = l)M(j 2 )l f \0 A2 \NN(IJ) > 

= 51 < j\j 2 'm\'m 2 \jm >< jlfmrrif\JM > 

j 

■\k\Y lfmf (k) < VM\\0 A2 \\NN(IJ) > (1) 
with the reduced matrix element defined as 

< VM\\0 A2 \\NN(IJ) >= F(k) <IJ -> VM > SF B(I, J) . (2) 

The atomic pp state is specified by isospin component I and total angular momentum J = 

0. 1, the latter values corresponding to the 1 S'o and 3 SDi states respectively. The two-meson 
state VM is specified by the intrinsic spin ji >2 , the total spin coupling j, the relative orbital 
angular momentum If = 1 and the relative momentum k. Eq. (0) includes a final state 
form factor F(k), the spin-flavor weight < IJ — > VM >sf and an initial state interaction 
coefficient B(I, J), containing the distortion in the protonium state J with isospin component 

1. Detailed expressions for these factors are summarized in Appendix 0. 

For the process V — > 7 (Fig. 2b), where the outgoing photon with energy k° is on-mass 
shell, we obtain, with the details shown in Appendix 0: 

T y ^ 7 = e ■ %i) Tr(Q Vv ) j^-^j ( 3 ) 

where e and S(mi), with projection mi, are the polarization vectors of 7 and V, respectively. 
The flavor dependence of the transition is contained in the factor Tr(Q(py), where Q is the 
quark charge matrix and ipv the QQ flavor wave function of vector meson V. In setting up 
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the two-step process NN(IJ) — > V M — > 7M we use time-ordered perturbation theory with 
the resulting matrix element ||T6f 



Tnn(IJ)-*VM-*jM — X! ^V->7 — 2 T NN (jj^ VM (4) 

mi m v s 

where the relativistic propagator for the intermediate vector meson in a zero width approx- 
imation is included. We resort to a relativistic prescription of the vector meson, since, with 
the kinematical constraint y/s = 0, V has to be treated as a virtual particle, which is severely 
off its mass-shell. Accordingly, an additional factor 2my, with the vector meson mass my, 
has to be included to obtain the proper normalization. From redefining 

Illy 6 

we generate the standard VDM expression of 

Tnn{ij)^vm^m — T N ft (I j)_> VM e ■ S(mi) A Vl . (6) 

mi 

The VDM amplitude Ay 1 , derived in the quark model, is: 



A Vl = V2 Tr(Qip v )^J-j^-y , (7) 

which in the limit m v Pd k° reduces to the well-known results of p] 

A n — e/fp — 0.055 and A^ = -A^ . (8) 

o 

The phase structure of Ay 7 , as determined by (pv, is consistent with the corresponding 
definitions for the strong transition matrix element. 

In the radiative annihilation amplitude, the coherent sum of amplitudes for V = p and 
u, arising from different isospin channels, has to be taken. This gives: 

Tnn(J)^M — X! ^ ' '-^'NN(IJ)—*VM-^'yM (9) 
V=p,u> 

where 5=1 for V ^ M and 5 = \/2 for V = M. The additional 6 accounts for the two 
possible contributions to the amplitude from an intermediate state with V = M, including 
a Bose-Einstein factor. For the decay width of NN — > 7X we write 



(10) 

M,e T ,m 2 V > 

Pf is the final state density and the sum is over the final state magnetic quantum numbers 
of meson X (1712) and of the photon (with transverse polarization ey). The corresponding 
branching ratio B is: 



where a statistical weight of the initial protonium state J with decay width T to t{J) is taken. 
With the details of the evaluation indicated in Appendix we finally obtain for the branch- 
ing ratios of pp — > 7X (X = 7r°, 77, 77'): 

U ; 41^(7=1) 7 W ' ; " 

•|B(0, 1) < 13 S£>x - P V > 5F 1) < 33 SD X -> c^vr > 5F | 2 . (12) 

Alternatively, .6(771-°) can be expressed in terms of the branching ratios .B(V7r ) for the 
strong transitions NN — > V7r° (Eq. ( |A12j ) of Appendix fS|): 

5(7^°) = T^sy^ (^(A ) + ^(^°) + low^V^C^W^ (is) 

with the interference phase f3j=\ determined by 

ife{g(0,l)*fi(l,l)} 
COS ^ =1 = |fi(0,l)B(l,l)| • (M) 

The same equations apply for X = rj and rf with 7r° being replaced by the respective meson. 

Here, a kinematical phase space factor / is introduced, which can be identified with those 

derived in specific models (Eqs. ( All|) and ( U10| )) or taken from phenomenology. Values 



for the branching ratios on the right hand side of Eq. ([13]) can either be taken directly 
from experiment or determined in the quark model considered in Appendix [A]. Magnitude 
and sign of the interference term, as determined by cosf3j=i, solely depends on initial state 
interaction for the spin-triplet NN state (J=l), as expressed by the coefficients £>(/, J = 1). 

Similarly, for the branching ratios of pp — > 7X (X = p ,u,j), now produced from the 
spin-singlet state (J=0) of protonium, we obtain: 



B(1P°) = jjj$ A % (1 B (P ^) + 2B (P°P°) + 2 ^cos{3 J= ^B{p^)B{p»p^ , (15) 
B{loo) = J^Ki ( B (P°^ + \ B b*>) + 2 -^cos(3 J= ^B{p^)B{^ , (16) 

and 

5(77) = jj0j A % { B + iB(u>fP) + ±B{uu) + +yB(pY)B(uJUj)+ 

+ *fcosP J=0 ^B(p°u) (^S(pPpP) + yB{uxj))} (17) 

with the interference determined as 

ite{g(0,0)*iB(l,0)} . 

cos ^° = \mwm\ ■ ( } 

Again, the sign and size of the interference cosf3j = o are fixed by the the initial state inter- 
action, here for protonium states with J=0. 

Eqs. (0), ((nj) and flTB] ) are analogous to those of Ref. 0; this is also true for Eq. (|T 



in the SU(3) flavor limit with B(p°p°) = B(uu). However, the essential and new feature 
of the present derivation is that the interference term is completely traced to the distortion 
in the initial protonium state. The possibility to link the interference terms cosf3j to the 
initial state interaction in protonium is based on the separability of the transition amplitude 
Tnn(ij)^vm- The sign and size of cos[3j (J=0,1) will have a direct physical interpretation, 
which will be discussed in the following chapter. 

We briefly comment on alternative model descriptions for the strong transition ampli- 
tudes iViV — > VM and its consequences for the interference terms in radiative pp decays. 
Competing quark model approaches in the description of NN annihilation into two mesons 
concern rearrangement diagrams as opposed to the planar diagram of the A2 prescription of 
Fig. 2a. In the rearrangement model a quark- ant iquark pair of the initial NN state is anni- 
hilated and the remaining quarks rearrange into two mesons. The quantum numbers of the 
annihilated quark-antiquark pair are either that of the vacuum ( 3 Pq- vertex, R2 model Ji"7|) 
or that of a gluon ( 3 Si-vertex, S2 model [18|,19|). In the R2 model, two ground state mesons 
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cannot be produced from an initial NN state in a relative S-wave; hence R2 is not applica- 
ble to the annihilation process considered here. The S2 model generates transition matrix 
elements for pp — ► VM, which are analogous to the ones of the A2 model of Eqs. ([I]) and 



(g), but with different absolute values for the spin-flavor weights < I J — ► VM >$f p8| , [19 
However, the relative signs of the matrix elements < I J — > pM > and < I J — > uoM > 
are identical to the ones of the A2 model, except in the case M = r\ where it is opposite. 
Therefore, results for branching ratios B{^M) of radiative decays expressed in terms of the 
branching ratios B(VM) are identical both in the A2 and the S2 approach, except for -8(7/7) 
where cosj3j=i changes sign. But, as will be shown later, the sign structure of cos[3j deduced 
in the framework of the A2 model is consistent with the one deduced from experiment. 

Possible deviations from the formalism presented here include contributions from virtual 
NA ± AN and AA states to the annihilation amplitudes as induced by initial state inter- 
action. The role of A state admixture and its effect on pp annihilation cross sections in 
the context of quark models was studied in Refs. [I3|,20j. Although contributions involving 



annihilation from iVA and A A states can be sizable POL the overall effect on the annihi 



lation cross section is strongly model dependent. In the case of the A2 model |T3L these 



contributions are found to be strong for NN D-wave coupling to channels with a virtual 
A in the S-wave, hence dominantly for the 13 SDi partial wave, where for isospin 1 = 
the tensor force induces strong mixing. However, for the radiative decay processes at rest 
considered here, the possible NA ± AA^ configurations only reside in the 33 SDi state (here 
33 SDi — > 7i°uj and 33 SDi — > rjp ). Due to the weak D-wave coupling in the 1=1 channel, 
NA configurations play a minor role and are neglected. 

Alternatively, the strong transition amplitudes NN — > VM can be derived in baryon 
exchange models Here however, the analysis is strongly influenced by the presence 

both of vector and tensor couplings of the vector mesons to the nucleon, by contributions of 
both N and A exchange (where the latter one contributes to the p°p° and 7r°p° channels) 
and by the addition of vertex form factors. The interplay between these additional model 
dependencies complicates an equivalent analysis. Due to simplicity we restrict the current 
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approach to a certain class of quark models, although deviations from the analysis given 
below when applying baryon exchange models cannot be excluded. 



III. PRESENTATION OF RESULTS 



In Sec. [Ill A] we discuss the direct application of the quark model approach to the radia- 



tive NN annihilation process. In Sec. |I I IB] we focus specifically on the isospin interference 
effects occuring in radiative transitions. We show that the the interference term is solely 
determined by the isospin dependent NN interaction, and give theoretical predictions for 
the phase cos/3j in various NN interaction models. Sign and size of cosflj can be interpreted 
by the dominance of either the pp or the nn component of the protonium wave function in 
the annihilation region. Furthermore we show that extraction of the interference term from 
experimental data is greatly affected by the choice of the kinematical form factor. Finally we 
comment on the applicability of the vector dominance approach to the pp — > 7^ transition. 



A. Branching ratios of radiative protonium annihilation 

In a first step we directly evaluate the expression for B{-K Qr y) and B(Xj), X = n, u, rf, p 
and 7, as given by Eq. fljjD , fljjD - fll7D and Eq. (|A12|) . To reduce the model dependencies 
we choose a simplified phenomenological approach as advocated in studies for two-meson 
branching ratios in NN annihilation ||21|| . The initial state interaction coefficients 13(1, J) 
are related to the probability for a protonium state with spin J and isospin I, with the 
normalization condition \B(0, J)\ 2 + |£>(1, J)| 2 = 1. The total decay width of state J is given 
by Ttot(J) with the separation into isospin contributions as T to t(J) = ^o(J) + Ti(J). We 
identify the ratio of isospin probabilities \B(0, J)\ 2 /\B(1, J)\ 2 with that of partial annihilation 
widths T (J)/Ti(J). For our calculations we adopt the isospin probabilities deduced from 
protonium states obtained with the Kohno-Weise NN potential |22] , where pp — nn isospin 
mixing and tensor coupling in the the 3 SD\ state are fully included |23|]. The resulting 
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values for B(I, J) are shown in Table |. The kinematical form factor f(j, X) is taken of the 
form ^| 

/(Tj X) = k- exp^-A^s-ml} (19) 

where k is the final state cm. momentum with total energy a/s. The constant A = 
1.2 GeV~ l is obtained from a phenomeno logical fit to the momentum dependence of var- 
ious multipion final states in pp annihilation |24]. Results for the branching ratios in this 
simple model ansatz are given in Table |TJ. For the decay modes 777 and 7/7 we use a pseu- 
doscalar mixing angle of O p = —17.3° ||25|| . The model contains a free strength parameter, 
corresponding to the strong annihilation into two mesons in the two-step process. Since we 
compare the relative strengths of the branching ratios, we choose to normalize the entry for 
-6(771"°) to the experimental number. The A2 quark model prediction for the hierarchy of 
branching ratios is consistent with experiment. In particular, the relative strength of tran- 
sitions from the spin-singlet ( 1 5'o) and triplet ( 3 SDi) NN states is well understood. The 
results of Table [Djgive a first hint, that the VDM approach is a reliable tool in analysing the 
radiative decays of protonium. Furthermore, all considered branching ratios are consistent 
with minimal kinematical and dynamical assumptions. We stress that the good quality of 
the theoretical fit to the experimental data of Table || should not be overemphasized given 
the simple phenomenological approach where initial state interaction is introduced in an 
averaged fashion. Although the A2 model provides a reasonable account of NN annihila- 
tion data, discrepancies remain in certain two-meson channels |I,2l| . In particular, observed 



two-meson annihilation branching ratios can show strong deviations from simple statistical 
or flavor symmetry estimates (dynamical selction rules), which in their full completeness 
cannot be described by existing models. Furthermore, theoretical predictions for two-meson 
branching ratios can be strongly influenced by initial state NN interaction (see for example 
Ref. [PB|), as in the case of radiative decays, but also by the possible presence of final state 
meson-meson scattering [|T5],p7]]. Given these limitations in the understanding of two-meson 



annihilation phenomena we will in turn dominantly focus on the determination of the inter- 

11 



ference term present in radiative pp decays. Here NN annihilation model dependencies are 
avoided by resorting to the experimentally measured two-meson branching ratios. 



B. Isospin interference and initial state interaction 

In a second step we focus on the determination and interpretation of the isospin inter- 



ference terms cosf3j (J=0,1) given by Eqs. (|lj) and fll8|), which in turn are related to the 
NN initial state interaction via the coefficients B(I, J) in Eq. (|A8|) . 

A full treatment of protonium states must include both Coulomb and the shorter ranged 
strong interaction, where the coupling of pp and nn configurations is included. The isospin 
content of the corresponding protonium wave function \1/ depends on r; for large distances ^ 
approaches a pure pp configuration, i.e., *ff(I = 0) = — 1). As r decreases below 2 fm, \l/ 
starts to rotate towards an isospin eigenstate, i.e. \& takes the isospin of the most attractive 
potential in the short distance region. The AW annihilation process under consideration 
here is most sensitive to the behaviour of \I> for r < 1 fm, where the strong spin- and isospin 
dependence of the NN interaction may cause either the I = or the 1=1 component 
to dominate. The consequences of the spin-isospin structure for energy shifts and widths 
of low lying protonium states have been discussed in Refs. [p3[ , p8l , ^9H ■ The sensitivity of 



pp — nn mixing in protonium states to changes in the meson-exchange contributions to the 
NN interaction was explored in [f2~l|1 . 

Let us first discuss the physical interpretation of the interference terms cos(3j. For a 
protonium state described by a pure pp wave function, the isospin dependent initial state 
interaction coefficients are related, B{1 = 0, J) = E{1 = 1, J). Similarly, for a protonium 
state given by a pure nn wave function in the annihilation region, that is ^>(I = 0) = 



= 1), B(I = 0, J) = -B(I = 1, J). Together with Eqs. flTJ) and (|18|), we obtain for 
the interference terms 

„ +1 for pure pp 
cos(3j = { . (20) 

-1 for pure nn 
12 



Therefore, a dominant pp component in the protonium wave function leads to constructive 
interference in radiative annihilation, with cos(3j = 1. Destructive interference reflects the 
dominance of the nn component in the annihilation region of the protonium state. Given 
this direct physical interpretation of the interference terms, radiative annihilation serves as 
a indicator for the isospin dependence of the NN protonium wave functions. 

For quantitative predictions of the interference terms and for comparison, we resort to 
protonium wave functions calculated p3f with three different potential models of the NN 
interaction, that by Kohno and Weise p2f (KW) and the two versions of the Dover-Richard 
30| , |29|| (DR1 and DR2) potentials. The calculation of Ref. [23] takes fully account of the 



neutron-proton mass difference, tensor coupling and isospin mixing induced by the Coulomb 
interaction. 

Results for the interference terms cos(3j as deduced from the three different potential 
models are given in Table |TJ. The value of the range parameter cIaz in the initial state 
form factor entering in Eq. ( |A8| ) is adjusted to the range of the annihilation potential of the 
respective models. With the choice of d,A2 — 0.12 fm 2 (KW and DR2) and g?a2 = 0.03 fm 2 
(DR1) the calculated ratios of isospin probabilities \B(0, J)\ 2 /\B(1, J)\ 2 are close to those of 
partial annihilation widths Tq(J) /T\(J) calculated in Ref. PBf . All three potential models 



consistently predict constructive interference for radiative annihilation from the atomic 1 Sq 
state, indicating a dominant pp component. For radiative annihilation from the spin triplet 
state 3 Si predictions range from nearly vanishing (DR1) to a sizable destructive interference, 
where latter effect can be traced to a dominant short ranged nn component in the protonium 



state. In Table |TTTj we also indicate predictions for the interference term cosfii, where the 
D-wave admixture in the 3 SDi state has been included. The results are obtained for the 
specified values of g^ 2 with an additional choice of hadron size parameters (that is R 2 N / R 2 M = 
0.6 or < r 2 >]{ 2 / < r 2 > 1 m = 1.2) entering in the expression of Eq. (|A8| ) . The inclusion 
of D-wave admixture in the initial state interaction coefficients 13(1, J = 1) as outlined in 
Appendix ^ is a particular feature of the A2 quark model. Hence, predictions for cosj3\ 
with the 3 Di component of the atomic 3 SD\ state included are strongly model dependent 
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and should not be overestimated. Generally, inclusion of the D-wave component in the form 
dictated by the quark model tends to increase the values of the interference terms. 

We also investigated the sensitivity of the interference term cos/3 j on the range of the 
initial state form factor, expressed by the coefficient d,A2- Although the absolute values for 
the initial state interaction coefficients 13(1, J) sensitively depend on the specific value for 
Aa2i variation of dA2 by up to 50 % has little influence on sign and also on size of cosfij. 
Thus, predictions for the interference terms cosfij in all three potential models considered, 
are fairly independent on the specific annihilation range of the NN initial state. 

The models used for describing the NN initial state interaction in protonium are char- 
acterized by a state independent, complex optical potential due to annihilation. Potentials 
of this type reproduce the low-energy pp cross sections and protonium observables, such as 
energy shifts and widths, fairly well. A more advanced fit to NN scattering data, in 
particular to the analysing powers for elastic and charge-exchange scattering, requires the 
introduction of an explicit state and energy dependence in the phenomenological short range 
part of the NN interaction. At present, latter NN potential |J]| was not applied to the 
protonium system; hence the model predictions of Table [ITJ should be regarded as a first 
estimate for the pp — nn mixing mechanism in the NN annihilation region. 



C. Isospin interference from data 

The VDM approach allows to relate the branching ratios of radiative annihilation modes 
to branching ratios with final states containing one or two vector mesons. Using these 
measured branching ratios in Eqs. ( |13D and ([15]) - (|17]) we can extract the interference 



terms cosftj directly from experiment. However, conclusions on the sign and size of the 
interference terms strongly depend on the choice of the kinematical form factor f(X\, X2), 
X\ and X2 = 7 or meson, entering in the different expressions. A first analysis || for 
determining the interference terms from data was performed by the Crystal Barrel Group, 
assuming a form factor of 
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where k is the final state cm. momentum and the interaction range is chosen as R = 1.0 fm. 
This form factor is appropriate for small momenta k, taking into account the centrifugal 
barrier effects near threshold. However, for radiative decays, with high relative momenta 
in the final state, the choice of Eq. QTS|) is more appropriate, it contains an exponential 
which restricts the importance of each decay channel to the energy region near threshold. 
This can be regarded as a manifestation of multichannel unitarity, that is the contribution 
of a given decay channel cannot grow linearly with k (as in the form of Eq. fl2lD ), since 
other channels open up and compete for the available flux, subject to the unitarity limit. 
Also, the latter form factor is given a sound phenomenological basis in NN annihilation 
analyses, for a more detailed discussion see for example Ref. ||33|| . Extracted values for the 



interference terms cos(3j with different J and different prescriptions for the kinematical form 
factor are given in Table [IV]. We also include there a third choice for the kinematical form 
factor (Eq. ( |A11| )), as deduced from the A2 quark model description of the NN annihilation 
process. Although finite size effects of the hadrons are included here, through the harmonic 
oscillator ansatz for the hadron wave functions the form factor is again useful for low relative 



momenta k. For the results of Table [TV] we use the measured branching ratios of pp — > ti p 



|34j1 , 7r°u, rjuj, ujuj, rj'uj H35] , rjp f36j or |j37H , HH an d rfp |33|. Values for cos(3j using the 



phase space factor of Eq. (|2lD are directly taken from the original analysis of Ref. ||. Error 
estimates for the other entries in Table [Iy| assume statistical independence of the measured 
branching ratios. For the radiative decay channel 7/7 only an upper limit for cos(3\ can be 
given. 

For all three choices of the kinematical form factor, the extracted values of cosf3j are 
consistent with the VDM assumption as they correspond to physical values. However, as 
evident from Table [TV], conclusions on sign and size of the interference strongly depend on 
the form of the kinematical phase space factor. For the preferred choice, i.e. Eq. (jT9|), we 
deduce destructive interference for radiative annihilation from the 3 SDi state, while for the 
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1 S'o state the corresponding isospin amplitudes interfere constructively. This is in contrast 
to the original analysis of Ref. ||, where the interference term is determined to be almost 
maximally destructive for all channels considered. Given the large uncertainties for cosftj 
using the preferred form factor, the values deduced from data are at least qualitatively 
consistent with the theoretical predictions of Table [TTT| , indicating a dominant pp component 
for the 1 S'o and a sizable nn component for the 3 SDi protonium wave function. As discussed 



in Sec. [I I fB], precise values for cosfij are rather sensitive on the isospin decomposition of the 
protonium wave function in the annihilation region. However, the current uncertainties in 
the experimental data should be very much improved to allow a more quantitative analysis 
of the isospin dependence of the NN interaction. 

D. Vector dominance model and the pp — > j(f> transition 

Measurements on nucleon-antinucleon annihilation reactions into channels containing (f> 
mesons indicate apparent violations of the Okubo-Zweig-Iizuka (OZI) rule f33f . According to 



the OZI rule, <fi can only be produced through its non-strange quark- ant iquark component, 
hence <fi production should vanish for an ideally mixed vector meson nonet. Defining the 
deviation from the ideal mixing angle 6 = 35, 3° by a = 9 — 6 and asssuming the validity 
of the OZI rule, one obtains the theoretically expected ratio of branching ratios [|TJ: 

R(X) = B(NN -> <f>X) /B(NN -> ujX) = tan 2 a w 0.001 - 0.003 (22) 



where X represents a non-strange meson or a photon. Recent experiments [|33| have provided 
data on the <p/uj ratios which are generally larger than the standard estimate of Eq. (p2|). 
The most notable case is the 07 channel for pp annihilation in liquid hydrogen |10] , where 
data show a dramatic violation of the OZI rule of up to two orders of magnitude, that is 
R(X = 7) w 0.3. Substantial OZI rule violations in the reactions pp — > Xcf) can possibly 
be linked to the presence of strange quark components in the nucleon [ft9p0| . However, 
apparent OZI rule violations can also be generated by conventional second order processes, 



even if the first order term corresponds to a disconnected quark graph ||T].f|2 . 
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In Refs. |42| , |TT| the apparently large value for the branching ratio -8(70) for pp an- 
nihilation in liquid hydrogen is explained within the framework of the VDM. Using the 
experimental rates of B(p<p) = (3.4 ± 1.0) x 1(T 4 and B(ucf>) = (5.3 ± 2.2) x 10" 4 |3j as 
inputs, the branching ratio -8(70) is given in the VDM by: 

J{ UJV ) 

Since the (put and (pp channels also violate the OZI rule estimate, R(X = p) m R(X = uj) ~ 
10" 2 ||33|| , the standard uj— <p mixing cannot be the dominant mechanism for the production of 



the (puj and <pp channels and the formalism developed in Sec. |I| cannot be used to determine 
the phase structure of the interference term cos(3q for .8(70). Consequently, the interference 
term cos(5q extracted in the 70; reaction is not necessarily consistent with that of the 70 
decay channel. For maximal constructive interference (cos/3 = 1) one obtains an upper limit 
for -6(70) in the VDM calculation of: 

5(70) = 2.7 x 1(T 5 for f = k 3 [11] 

5(70) = 1.5 x 10" 6 for / given by Eq. (H) (24) 

This is to be compared with the experimental result -8(70) = (2.0 ± 0.4) x 10~ 5 [|10|| . The 
possibility to explain the experimental value of -8(70) in VDM depends again strongly on 
the choice of the kinematical form factor. In Ref. |ll| the form / = k 3 is used, appropriate 
for relative momenta k near threshold, resulting in an upper limit which lies slightly above 
the observed rate for 5(70). With the choice of Eq. ([T9| ) the upper value underestimates 



the experimental number by an order of magnitude. When we extract the interference 
terms cosj3j from the conventional radiative decay modes with the choice / = k 3 , we ob- 
tain: cosfii = —1.32 for 7T7, cosfii = —0.94 for 177 and cosj3 = —.90 for uj. Hence, a 
near threshold prescription for the kinematical form factor in the VDM leads to maximal 
destructive interference for all channels considered, exceeding even the physical limit in the 
case of 717. This would indicate a nearly pure nn component in the annihilation range of the 
protonium wave functions for both the J=0 and 1 states. These results are in strong conflict 
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with the theoretical expectations for cosfij reported in Sec. |I 1B|, where at least qualitative 



consistency is achieved with the kinematical form factor of Eq. (O) . 



Recent experimental results [0] for the reaction cross section pp — > 00 exceed the simple 
OZI rule estimate by about two orders of magnitude. Therefore, in the context of VDM 
an additional sizable contribution to the branching ratio -6(70) might arise, although off- 
shell, from the 00 intermediate state. With an estimated cross section of pp — > uouj of 
about 5 mb in the energy range of the 00 production experiment, the ratio of cross sections 
is given as a^ja wul ~ 3.5 //6/0.5 mb Given the measured branching ratios of ujuj 



| 35| and tocj) we can simply estimate the ratio of strong transition matrix elements for 



annihilation into 00 and uocf) from protonium of y_B(00)/S(a;0) « 0.43. For this simple 
order of magnitude estimate we assume that a^/ a UUJ is partial wave independent and phase 
space corrections are neglected. With the VDM amplitude = ^A P1 we obtain an 
upper limit of -6(70) ~ 2.3 x 10 -6 with f given by Eq. (|T9|), where the contribution of the 
00 intermediate state is now included. Excluding an even further dramatic enhancement 
of the 00 channel for NN S-wave annihilation, inclusion of the 00 intermediate state does 
not alter the conclusions drawn from the results of Eq. (0). Hence, the large observed 
branching ratio for 70 remains unexplained in the framework of VDM. 



IV. SUMMARY AND CONCLUSIONS 

We have performed a detailed analysis of radiative pp annihilation in the framework 
of a two-step process, that is pp annihilates into two-meson channels containing a vector 
meson which is subsequently converted into a photon via the VDM. Both processes are 
consistently formulated in the quark model, which allows to uniquely identify the source 
of the isospin interference present in radiative transitions. Based on the separability of the 
transition amplitude NN — > VM, sign and size of the interference terms can be linked to 
the dominance of either the pp or the nn component of the Is protonium wave function in 
the annihilation region, hence constitutes a direct test of the isospin dependence of the NN 
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interaction. 

In a first step we directly applied the quark model in a simplified phenomenological 
approach to the radiative NN annihilation process. Model predictions are consistent with 
data and confirm the usefulness of VDM in the analysis of radiative transitions. In a 
second step we discussed sign and size of the interference term as expressed by cos(3j (J = 
0,1). Direct predictions of cosfij, as calculated for different potential models of the NN 
interaction, are qualitatively consistent, in that a sizable constructive interference is deduced 
for radiative annihilation from the atomic x Sq state, while for the 3 Si state the interference 
term is vanishing or destructive. These predictions should be tested with more realistic 
parameterizations of the NN interaction ||31|| . Extraction of the interference effect from 
data is greatly influenced by the choice of the kinematical form factor associated with the 
transition. Values of cos/3j determined for the preferred form of Eq. (|19|) are qualitatively 
consistent with our theoretical study; however, a more quantitative analysis is restricted by 
the present uncertainties in the experimental data. Within the consistent approach emerging 
from the analysis of non-strange radiative decay modes of protonium, an explanation of the 
measured branching ratio for the OZI suppressed reaction pp — > 7$ cannot be achieved. New 
mechanisms, linked to the strangeness content in the nucleon, may possibly be responsible 
for the dramatic violation of the OZI rule in the 7$ final state. 
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APPENDIX A: NUCLEON-ANTINUCLEON ANNIHILATION INTO TWO 

MESONS IN THE QUARK MODEL 



In describing the annihilation process of iViV — > VM where V = p, u and M = 
7T°, 7], p, uj and rj we use the A2-model of Fig. 2a. Detailed definitions and derivation 
of this particular quark model are found in Refs. [P],|13[. The initial state A^A" quantum 
numbers are defined by i = ILSJM (I is the isospin, L is the orbital angular momentum, 
S is the spin and J is the total angular momentum with projection M). For the final two 
meson state VM we specify the angular momentum quantum numbers, with j 12 indicating 
the spin of mesons 1 and 2, j the total spin coupling and If the relative orbital angular mo- 
mentum. For the transitions of interest the quantum numbers are restricted to L=0 and 2, 
corresponding to pp annihilation at rest in liquid hydrogen, ji = 1, representing the vector 
meson, and If = 1, given by parity conservation. Taking plane waves for the initial and 
final state wave functions with relative momenta p and k, respectively, the transition matrix 
element is given in a partial wave basis as: 

Tnn^vm = < V( 3l )M( 32 )lf\0 A2 \NN{i) > 

= < jij2fni^2\jfn >< jlfmm,f\JM > 
j 

■\k\Y lfmf Ck)Y^ s M \p) < VM{j u j 2 ,j,l f )\\0 A2 \\NN(i) > . (Al) 

The reduced matrix element of the two-meson transition is given in the A2 model as: 

< VM\\0 A2 \\NN{i) >= F Wf p L exp(-d A2 (3/Ak 2 + l/3p 2 ) <i^VM> SF . (A2) 

The factor i 7 L,j / is a positive geometrical constant depending on the size parameters of the 
hadrons for given orbital angular momenta L and If. The exponentials arise from the overlap 
of harmonic oscillator wave functions used for the hadrons with the coefficient d A2 depending 
on the size parameters and Rm of the nucleon and meson: 

dA2 = (A3) 
6K N + ZK M 
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The matrix elements < i — > VM >sf are the spin-flavor weights of the different transitions 
listed in Table [V]. Note that with the flavor part of the vector mesons defined as 

p° = —j={uu — dd), uj = —=(uu + dd) (A4) 
V 2 v 2 

the matrix elements < 2 — > pM > and < i — ► cjM > have same sign. For the tensor force 
coupled channel 3 SDi the spin-flavor matrix elements are simply related by a proportionality 
factor, dependent on the isospin channel, but independent of the VM combination, that is: 



Fl-2,i,,i < 2,+vi Di -> VM > SF = C(I) • F^tt-i < 2 ' +1 ' S s i -> VM >sf , 

CO = (,!?' -*£) ' • (AS) 



In coordinate space the protonium wave function, including tensor coupling and isospin 
mixing, is written as: 

%p(J, 5) = E ^LsAr)Yj! s M (r) . (A6) 
Inserting this wave function into the expression for the transition matrix element results in: 

T N n(ij)-+vm = < ji32m 1 m 2 \jm >< jl f mm f \JM > 
j 

■\k\Y lfmf (k)F(k) <i^VM > SF B(I, J) , 
F(k) = exp(-d A2 3/4k 2 ) . (A7) 

The distortion due to initial state interaction is contained in the coefficient £>(/, J), which is 
simply the overlap of the isospin decomposed protonium wave function with the effective ini- 
tial form factor arising in the transition. By taking the Fourier transform of the initial state 
form factor contained in Eq. (|A2|) , these coefficients for the Is atomic states of protonium 
are defined as: 

B(I, J = 0) = F L=0 , lf=1 (2d A2 /3)- 3/2 jH drr 2 exp(-3r 2 /(4d A2 ))^ mo (r) for ^ , 

{POO 
(2d A2 /3)- 3/2 J Q drr 2 exp(-3r 2 /(4d A2 ))^ mi (r)- 

-C{I) (2d A2 /3)- 7/2 drr 4 exp(-3r 2 /(4^ 2 ))^ /211 (r)} for 3 SD 1 . (A8) 
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The partial decay width for the annihilation of a protonium state with total angular 
momentum J into two mesons VM is given by 

T pP ^ VM {I, J) — 27T V M k f dk \ t nn(ij)^vm\ 2 (A9) 



where E is the total energy and Ey jM = y rriy^ M + k 2 the energy of the respective outgo- 
ing meson with \k\ fixed by energy conservation. With the explicit form of the transition 
amplitude of Eq. (|A7|), the partial decay width is written as: 



T PP ^ VM (I, J) = f(V, M)<i^VM > 2 SF \B(I, J) | 2 (A10) 

with the kinematical phase space factor defined by: 

f(V, M) = 2n^^k 3 exp (-3/2d A2 k 2 ) . (All) 

Taking an admixture of initial states given by their statistical weight, the branching ratio 
of S-wave pp annihilation into the two meson final state VM is given by: 

B(VM) = B( P p - VM) = £ ^ + l)r^vM(/,J) (A12) 

j 41 tot[J) 
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APPENDIX B: VECTOR MESON - PHOTON CONVERSION IN THE QUARK 

MODEL 



The transition V — > 7 (Fig. 2b), where V = p or u, can be formulated in the quark 
model, and related to the physical process of V — > e + e~ . An explicit derivation of the latter 
process can be found in Ref. [JEJ. We just quote the main results necessary for the discussion 
of the radiative decays of protonium. 

The QQ'j interaction is defined by the Hamiltonian 



H I = e d^ m {x)AJx) (Bl) 



with the quark current 



= ff(5)Qyg(5) (B2) 



where q(x) is the quark field and A^(x) the electromagnetic field given in a free field ex- 
pansion. For emission of a photon with momentum k, energy k° and polarization e M from a 
vector meson with momentum py we obtain: 

<-y(k,eJ\Hr\V(fr) >=5(£-pV)Tv^ 7 (B3) 

with 

= -^jyr^ < < 0|j^(2 = 0)|V > . (B4) 

For the conversion of a vector meson V into a real photon only the spatial part of the 
current matrix element contributes. Using standard techniques for the evaluation of the 
current matrix element we obtain 

e\/ 6 -» 

Ty ^ = WW s Tr{Qvv) ^ = 0) (B5) 

with the quark charge matrix Q and the polarization S of the vector meson. The QQ flavor 
wave function (py is consistently defined as in Eq. ( [A4jj of Appendix ^ and contributes to 
the transition amplitude: 
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Tr(Q<p v )={ ^ H . (B6) 
573 for ^ 

The spatial part of the QQ wave function at the origin ip(r = 0) is given within the harmonic 
oscillator description as |?/>(0)| 2 = (ttR\, i )~ 3 ^ 2 , where the oscillator parameter Rm is related 



to to the rms-radius as < r 2 > 1 / 2 = ^J3/8R M . 

Extending the outlined formalism to the physical decay process V — > e + e~ the decay 



width is given as [45| 



IQna 2 r ... ^ . 1 2 



IW e - = — ^ {Tr(Q<^)} z |</>(0)| 2 (B7) 

with a = e 2 / (47r) and mass my of the vector meson. Latter result can be compared to the 
one obtained in the vector dominance approach resulting for example in 

47T a 2 m p 

T p o^ e+e - = — — — (B8) 
with the decay constant f p . Hence we can identify 

\m\ 2 = 7rf 2 (B9) 

which with the experimental result of r p o^ e + e - = 6.77 yields f p = 5.04 or equivalently 
Rm = 3.9 GeV -1 , very close to the preferred value obtained in the analysis of strong decays 
of mesons. Hence, the matrix element for the conversion of a vector meson into a real photon 
is alternatively written as: 

Ty^ = e- S Tr(Q<p v ) { ^Cf p ' ^ 
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APPENDIX C: MATRIX ELEMENTS AND DECAY WIDTH IN RADIATIVE 

ANNIHILATION 



In the following we present details for the evaluation of the matrix element of Eq. (6), 
which is explicitly written as: 

T NR{IJ) ^ VM ^ lM = J2 < 3i32mim 2 \jm >< jl f mm f \JM > \k\Y lfmf {k) 



mi 



■ < VM\\0 A2 \\NN{IJ) > e ■ S{m x )A v ^ (CI) 

— * 

where If = 1 and j = 1, for the processes considered. The relative final state momentum k 
and the photon polarization e are written in a spherical basis as: 

\k\Y lmf (k) = ]j^ k m f and e ■ %i) = e mi (C2) 

which together with Eq. (|(J1Q leads to the result: 



(C3) 



T N n(ij)^vm^m = VIf^7 < VM\\0 A2 \\NN(IJ) > ^ 
(e x k) M for j 2 = 0, J = 1 (M = vr , V ) 
^(e x £)_ m2 for j 2 = 1, J = (M = p°, w) 

Consequently, for the process iViV — > V\V 2 —>■ 7172 the transition matrix element is deter- 
mined as: 

T(NN(IJ) -> ViVa -> 7172) = E e m2 {2)Av 21 T{NN -> -> 71^2) 

m 2 

' A Vll A Vll < V\V 2 \\0 A2 \\N N {I J) > -^= (e(l) x fc) • e(2) (C4) 



where e(i) refer to the polarization of the photon i. 

The derivation of the decay widths for the radiative transitions is examplified here for 
the process NN — > 77?°. The corresponding matrix element is obtained by a coherent sum 
of intermediate vector meson states p and u> as: 

x fe)Af {^A < p 7r°\\O A2 \\ 13 SD 1 > +A U>1 < L07r°\\O A2 \\ 33 SD 1 >} . (C5) 
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The decay width for iViV — > 771" is then: 



T NR ^o = 2np f ^ ^TT |T(iV ^ (J) ~* 77T ° )|2 (C6) 



Pf = ^ldk, 1(7) 



with the final state density 

— # 

I A; I = k, and the sum is over the two transverse photon polarizations ex and the total 
projection M of the NN protonium with total angular momentum J. Using 

£ / dk\(ex k)\ 2 = 8irk 2 (C8) 

e T ,M J 

together with the expression for the reduced matrix element in Eq. (2), we finally obtain: 

r ArJ v^ 77r o = 

= /( 7 , 7r°)A 2 n \ < 13 SD l - pV > SF B(0, 1) + ^ < 33 SL>i - cuvr > 5F 1)| 2 (C9) 
with the kinematical phase space factor defined in analogy to Eq. ( |A11| ) as: 

/( T) M) = 2n^^exp (-3/2d A2 k 2 ) . (CIO) 
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FIG. 1. Two-step process NN -> MV -> M7 with F = p°, a; and M = ir°, rj, p°, uj, 17', 7 for 
radiative protonium annihilation. 

FIG. 2. Quark line diagrams corresponding to NN annihilation into two mesons (a) and vector 
meson - photon conversion (b). 
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TABLE I. Isospin probabilities \B(I, J)\ 2 - Values are deduced from calculation of partial anni- 
hilation widths of Is protonium with Kohno-Weise potential [23]. 



State |B(0, J)\ 2 \B(1, J)\ 2 T tot (J) [keV] 



% (J=o) 


0.60 0.40 


1.26 


3 5Di (J=l) 


0.53 0.47 


0.98 


TABLE II. Results 


for branching ratios B for pp — > 7 V with X = 7r°, r/, p°, 


7/ and 7 in the 


simple model estimate. The entry for i?(7r°7) is normalized to the experimental value. Data are 


taken from Ref. H. 






Channel 


5 x 10 6 (model) 


5 x 10 6 (exp.) 


3 S£>1 -> 7T°7 


44 


44 ±4 


3 ST>i -> rry 


14 


9.3 ± 1.4 


1 S Q -> ivy 


68 


68 + 19 


3 SD 1 -> 7/7 


8.3 


< 12 


-> 77 


0.14 


< 0.63 


-> P7 


50 





TABLE III. Isospin interference terms cos(5j as calculated with the Is protonium wave functions 
of the KW, DR1 and DR2 potential models. Values in brackets denote the results where the 3 D\ 
component of the atomic 3 SDi state is included with admixture fixed by A2 model. 





KW 


DR1 


DR2 


cos/?o 


+1.00 


+0.83 


+0.63 


cos(3\ 


-0.90 (-0.76) 


+0.10 (+0.36) 


-0.41 (+ 0.53) 
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TABLE IV. Isospin interference terms cosf3j as deduced from data. The label HQ refers to the 
kinematical form factor of von Hippel and Quigg as defined in Eq. (^T|). Similarly, labels VAN and 
A2 refer to the form factor prescription of Vandermeulen and of the A2 quark model, as defined 



in Eqs. ( |19|) and ( All ), respectively. The analysis for cosjij (HQ) is directly taken from Ref. [||. 
The first line of the analysis for 777 is done for B(r]p°) = (0.53 ± 0.14) x 10~ 2 ||], the second line 
for B(r]p°) = (0.33 ± 0.09) x 10~ 2 @. 

Channel cos(3j (HQ) cosf3j (VAN) cos/3 j (A2) 

3 SD 1 -> vr°7 ( J = 1) -0.75 ±0.11 -0.10 ±0.28 ±1.00 ±0.38 

3 SD 1 -> 777 ( J = 1) -0.78 ±0.25 -0.47 ±0.72 -0.43 ±0.76 

-0.58 ±0.48 -0.17 ±0.90 -0.12 ±0.96 

1 S ^u>y (J = 0) -0.60 ±0.18 ±0.15 ±0.38 -0.21 ±0.28 
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SD 1 -> r/7 (J = 1) < -0-26 < 1.65 < -0.12 



TABLE V. Spin-flavor matrix elements < i —> VM > S f for the decay NN(L = 0) — ► KM in 
the A2 quark model. These are relative matrix elements obtained from Ref. Q. Here, 77,^ refers 
to the nonstrange flavor combination 77^ = (uu ± dd) / y/2. 

Decay channel < i —* KM >s_f 



n S - 




-\/243 


11 

Do - 


-/V 
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Fig. 1 
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